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Abstract 

We measure the branching fraction of the radiative B~ decay B{B~ (j)K~'y) = (3.46 ± 
0.57+o;37) X lO"*', and set an upper hmit on the radiative B° decay ^ (/)F°7) < 2.71 x 10"^ 

at 90% confidence level. We also measure the direct CP asymmetry of the B~ (j)K~j mode 
Acp = (—26.4 lb 14.3 lb 4.8)%. The uncertainties are statistical and systematic, respectively. These 
measurements are based on 207 fb~^ of data collected at the T(4S) resonance with the BABAR 
detector. 
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Measurements of the branching fractions and CP asymmetries of 6 — > 57 decays provide a 
sensitive probe of the Standard Model (SM). In the SM these decays are forbidden at tree level 
but allowed through electroweak penguin processes (Fig. They are therefore sensitive to the 
possible effects of new physics ^ in the form of new heavy particles contributing to the loop 
diagram. Additional contributions to the decay amplitudes could affect branching fractions and 
CP violation. Furthermore, the radiated photon is polarized due to the left-handed nature of the 
weak interaction. The polarization can be probed by measuring the time-dependent CP asymmetry, 
which is sensitive to interference between B^-B^ mixing and decay. Theoretical estimates in the 
SM [21 bound the mixing-induced CP asymmetry at about the 10% level. Here we focus on the 
time-integrated direct CP asymmetry, which is expected to be the same for charged and neutral B 
decays. 

Although exclusive 5 — > 57 decays present a theoretical challenge due to large non-perturbative 
QCD interactions, they are experimentally clean. There have already been results published for 
branching fraction and/or CP asymmetry measurements in several exclusive modes: B — > K*^ [3], 
B^ K^t:^^ I^, 5 — > ijOK'j and various B Kinr'y JBl modes. Here, we present a 
measurement of the branching fraction for the charged mode B~ (f)K~'y and the neutral mode 
B^ 4>K^^. ^ We also measure the direct CP asymmetry in the charged mode Acp = — 
N{B^)\/[N{B^) + N{B~^)], where the flavor of the B is determined by the charge of the kaon. The 
Belle Collaboration has previously measured the branching fractions for these modes, using 90fb~^ 
of BB data at the T(4S') resonance [7]. We describe the first BABMi measurements of these modes 
using a dataset that is more than twice as large. 

The data used in this analysis were collected with the BABMi detector at the PEP-II asymmetric 
e^e~ storage ring. This analysis is based on a data set of 207 fb~^ corresponding to 228 million BB 
pairs collected at the T(4S') resonance. The BABAR detector is described in detail elsewhere [5]. 
Important to this analysis are the tracking system composed of the silicon vertex tracker (SVT) 
and drift chamber (DCH), the detector of internally reflected Cherenkov light (DIRC), and the 
electromagnetic calorimeter (EMC). The SVT and DCH provide tracking and ionization energy loss 
(dE/dx) measurements for charged particles inside a 1.5 T magnetic field. The SVT is composed 
of five layers of double sided silicon strips and covers a polar angle range between 20.1° and 150.2°. 

^Throughout this paper, whenever a mode is given, the charge conjugate is also implied. 
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The DCH continues tracking outside the SVT volume. It consists of 40 layers of hexagonal cells 
filled with an 80:20 mixture of helium:isobutane. The DIRC is a ring imaging Cherenkov light 
detector. Below 700 MeV/c the tracking system provides most of the charged particle identification 
(PID) information, while the DIRC contributes more information at higher momenta. Photons are 
detected and their energy measured in the EMC, which is composed of 6580 thallium-doped Csl 
crystals. 

We select events with EMC clusters of energy 1.5 — 2.6 GeV in the e~^e~ rest frame (CM 
frame) that are not associated with any charged track. Photon candidates are required to have an 
energy distribution consistent with the shower shape of an electromagnetic interaction, and they are 
required to be well isolated (> 25 cm) from other calorimeter clusters. A veto is applied to photon 
candidates that can be combined with other neutral EMC clusters above a minimum threshold 
energy to form an invariant mass consistent with a vr'' (115 — 155 MeV/c^) or rj (470 — 620 MeV/c^). 
The threshold energy is 50 MeV for vr'^ and 250 MeV for rj. 

We select (j) candidates from pairs of oppositely charged kaon tracks, determined not to be pions 
based on a PID likelihood selection algorithm that uses dE/dx and Cherenkov light measurements. 
The same algorithm is used for the single in the charged mode. The tracks are fitted to a vertex 
using a Kalman decay chain fitter IS] , and are required to have a vertex probability greater than 
0.1%. We select candidates with masses within a 10 MeV/c^ window of the nominal (j) mass |10| . 
In the neutral mode, pairs of oppositely charged tracks are fitted to a common decay vertex and 
accepted as Kg candidates if the fit yields a probability greater tha 0.1%. The invariant mass of 
the pair is required to be within 10 MeV/c^ of the Kg mass. The flight length of is required to 
be greater than three times the uncertainty of that length. We require the combined (pK invariant 
mass to be less than 3.0 GeV/c^. In the neutral mode a veto is applied by removing candidates 
with a (pK invariant mass within 10 MeV/c^ of the mass. 

The K, (p, and 7 candidates are fitted to a common vertex and accepted as B candidates 
if the vertex probability is greater than 0.1%. To discriminate BB events against continuum 
background the ratio of Legendre moments L2/L0 is required to be less than 0.55. The Li are 
defined by Li = J2j IpJI I cos where the p* are the CM momenta of all particles not used in 
reconstructing the signal B candidate, and the angle Oj is between the particle's momentum and 
the thrust axis of the signal B. We require the cosine of the angle between the B candidate and 
the beamline, cos0^, to be in the range [-0.9,0.9] in the CM frame. We use two uncorrelated 
kinematic variables of the B candidate: the reconstructed mass rrircc and the missing mass m^[ss- 
The reconstructed mass is the B candidate invariant mass calculated from the reconstructed energy 
and momentum. This is required to be within 4.98 — 5.48 GeV/c^. The missing mass squared is 
"^miss ~ (PBeams " pg**^^ '^°'^*'*') ^ , where pBeams IS the four-momentum of the beams and pg**^^ 
is the four-momentum of the B (pK'j candidate after a mass constraint on the B is applied. We 
require the missing mass to be in the range 5.12 — 5.32 GeV/c^. 

To study event selection criteria we use simulated Monte Carlo (MC) events of signal, generic B 
decays, and e+e~ — > qq continuum background, where q = {u,d,s,c}. Events are generated using 
EVTGEN ^2 and the detector response simulated with GEANT4 |12j . Signal Monte Carlo is generated 
according to the inclusive 6 ^ S7 scheme presented in reference jJSl; using nib = 4.62 GeV/c^ 
for the effective b quark mass. Exclusive signal MC is derived from this by using only the part 
of the hadronic mass spectrum above the (pK threshold of 1.52 GeV/c^. Our selection criteria 
were chosen to optimize the figure of merit Ns / \/ Ns + Nb in the signal region, where Ns and Nb 
are the number of signal and background events, respectively, and the signal region is defined by 
5.05 < mRec < 5.4 GeV/ c2, 5.27 < m^iss < 5.29 GeV/ c^, | cos 6*151 < 0.8, and L2/L0 < 0.48. 
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After all selection criteria are applied the average candidate multiplicity in events with at least 
one candidate is approximately 1.01 and 1.07 in the neutral and charged modes respectively. If 
multiple B candidates are found in an event, we select the best one based on a formed from the 
value and uncertainty of the mass of the (j) candidate and, in the neutral mode, the Kg candidate. 
The remaining background comes from continuum combinatorics and from other B decay modes, 
which will be discussed later. 

We use an extended nia,xinium likGliliood. fit in four observablcs ~ ''^miss? ^rec? 

L2/L0, and 

cos^Ij - to extract the signal and combinatoric background yields. The likelihood C is defined in 
the following way: 

^ = nf [{NsVh + NBVh)] . (1) 

Ns and Nb are the fitted number of signal and background events. is the total number of 
events used in the fit. and are products of the signal and background probability density 
functions (PDFs) for each event i. In the charged mode, in order to fit the CP asymmetries of 
the signal and the background, the numbers of B^ and B~ events are determined separately: 
Nj = ^(1 + fAcp)nj, where j = S 01 B; f is the fiavor, defined as +1 for B~ and —1 for B^; nj 
and Aqp are the total yield and CP asymmetry of species j. In the neutral mode A'^- = nj. 

The signal PDFs for m-miss and parametrized as asymmetric, variable-width Gaussian 

functions: 



/(x) = exp 



(2) 



The parameters (Jl r aiid a^^R determine the core width and variation of the width on either side 
of a; = 0. The mmiss background PDF is an ARGUS function P^, with the endpoint calculated 
on an event-by-event basis from the beam energy. The mrec background PDF is modelled as a 

figgj-gg polynomial. The L2/LQ distribution is modelled using a binned PDF with eight bins, 
because there is no a priori model for this distribution. There are seven parameters in the PDF due 
to the condition that the bins sum to unity. The signal and background models both use this form. 
The cos 6*^ distribution is modelled as a 2^^*^ degree polynomial in both signal and background; for 
true B candidates it is expected to follow 1 — cos^ 9*^. 

We use a high statistics B^ K*^{-^ K^ti^)^ control sample to determine our signal shape 
parameters. Once determined, these signal parameters are fixed for the fit to — > (j)K^ data, 
while all background shape parameters are allowed to vary. We fit for the number of signal and 
background events, and in the charged mode the signal and background CP asymmetry as well. 
The same fit is used with signal MC to determine the efficiency of the previously described selection 
criteria. Corrections to the efficiency are discussed below. 

We apply several corrections to the fitted signal yield and efficiency before determining the 
branching fractions. Studies of simulated events show that our main sources of peaking backgrounds 
are nonresonant B KK~^K~^ events, and B (pKir^, B 4>Kri, where the tt^ or r/ decay fakes 
a high energy photon. We estimate the amount of -B — > KK~^K~'y contamination by fitting for the 
yield in (p mass sideband regions extending outside the signal region from 10 MeV/ to 30 MeV/ 
of the nominal (f) mass. By interpolating into the signal region, we find and correct for 0.03 it 1.5 and 
5.4 lb 4.2 events for the neutral and charged modes respectively. These contributions are subtracted 
from the event yield found in the fit. We also subtract the expected amount of B ^ (f)K*{-^ Ktt^) 
as determined by BABAR |1^: 0.27 neutral and 1.98 charged events. Because there have been no 
branching fraction measurements of -B — > (pKir^, B — > (pKrj, we assume that the branching fraction 
of these modes is no more than three times that oi B ^ tpK* . Therefore, we assign an uncertainty 
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Channel 


Yield 


Efficiency 


e(io-«) 


Acp 


B- (pK'-f 


85.0 ± 13.91^:^ 
8.0± 5.5l:J:^ 
< 16.0 


[21.9± 1.6(syst)]% 
[15.33 ±0.81(syst)]% 


3.46 ± 0.57^^:1]^ 
1.35 ± 0.92t°JJ 
< 2.71 


(-26.4 ± 14.3 ±4.8)% 



Table 1: Summary of the branching fractions and direct CP asymmetry. In B^ (pK^^ the 90% 
confidence level upper limit is also given. 



of 0.51 neutral and 2.86 charged events due to nonresonant B (j)K{7r^/r]) background. To correct 
for any fit bias, we generate 1000 pseudo-experiments using our maximum likelihood PDFs with 
separate components for BB and continuum, and embedding signal events from the full simulation. 
The background components are generated using shape parameters determined from the full MC 
simulation. We correct for a bias of +4.07 it 0.45 events in the charged mode, due to correlations 
among the observables in signal MC events that are not accounted for in the fit. In the neutral 
mode we find a bias of —0.06 it 0.20, and so we include 0.20 events in the systematic uncertainty 
of the yield. We correct for known efficiency differences between data and Monte Carlo in charged 
track, single photon, and Kg reconstruction. These corrections amount to 0.956 in the neutral 
mode and 0.975 in the charged mode. The corrected efficiencies are (15.3 it 0.81)% in the neutral 
mode and (21.9 it 1.6)% in the charged mode, where the uncertainties are systematic (discussed 
below) . 

The signal yields, efficiencies, branching fractions, and charged mode CP asymmetry are re- 
ported in Tabled We calculate the central value of the branching fractions as: 

where A^sig is the corrected number of signal candidates, = (228.3 it 2.5) x 10^ is the number 
of BB pairs recorded by BABAR, e is the corrected efficiency, and the {Kg — > tt^ vr") term is 
only used in the neutral mode. The partial branching fractions are given by Ref. pOl. We measure 
B{B^ 4>K^-i) = (1.35±0.92+°JJ) x 10"^ and B{B~ (^K-j) = (3.46±0.57l[]J^) x 10"^. In the 
charged mode we measure Acp = (—26.4+14.3+4.8)%. Fits to the missing mass and reconstructed 
mass distributions, projected into the signal region defined earlier, are shown in Figure |51 We use 
a set of 1000 pseudo-experiments in the neutral mode to determine the probability of obtaining a 
branching fraction less than or equal to our measured central value under the hypothesis that it is 
in fact the same as the charged mode. This was found to be 1.1%. 

For the neutral mode we compute the 90% confidence level upper limit on the branching fraction. 
We use a Bayesian approach with an a priori probability for the branching ratio which is flat in the 
physical region < < 1, and zero elsewhere. The value of the likelihood function is computed by 
fixing the signal yield to a desired value and fitting the other free parameters on the data sample. 
The function is then integrated numerically. We account for systematic uncertainties on the yield 
by convolving the likelihood function with the distribution of the errors, before computing the 
upper limit. To determine the upper limit of the signal yield we use a Gaussian PDF having a 
width equal to the systematic uncertainty of the yield. Similarly for the efficiency uncertainty we 
use a Gaussian PDF having a width equal to the systematic error. After also applying the yield 
corrections discussed previously we obtain B{B^ — > (pK^j) < 2.71 x 10~^. 

To determine the contribution from resonances decaying to (j)K we study the background- 
subtracted UH]) efficiency corrected (j)K invariant mass distribution, shown in Figure |31 Using the 
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Figure 2: Missing mass (left) and reconstructed mass (right) fit projections in the signal region 
for the charged mode (upper) and the neutral mode (bottom). The dotted curves show the fitted 
background contribution while the solid curves show the signal. 
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Figure 3: The background-subtracted and efficiency-corrected (j)K mass distributions (points with 
uncertainties) for the charged mode (left) and the neutral mode (right). The signal MC prediction 
for the mass spectrum, based on Ref. ^H], is shown as a histogram without uncertainties. 

charged mode, we find that no more than 50% of the spectrum in the 1.6 — 3.0 GeV/c^ range comes 
from the -fC2(1770) resonance. We use this to bound our model uncertainty, described below. 

We assign an uncertainty due to the fixed signal parameters in the fit. The parameters were 
obtained from the control sample, and therefore they have some statistical uncertainty. We varied 
these parameters within their uncertainties to determine the total uncertainty on the yields. We 
account for other systematic uncertainties in charged kaon tracking, kaon PID, K'^, cj), and photon 
selection efficiency. There are small uncertainties assigned with our L2/-L0 selection and the vr*^/ 
r/ veto. To account for uncertainty due to the assumption of a specific (t)K mass spectrum for 
simulated events, we determine what our efficiency would have been in the scenario that half of 
the spectrum comes from i^2(1770) resonant (pK production, while the other half comes from our 
signal MC model. We assign the relative efficiency difference as an uncertainty. Adding all of the 
uncertainties in quadrature, we find a total acceptance/efficiency uncertainty of 5.2% in the neutral 
mode and 7.1% in the charged mode. The contributions for each mode are summarized in Tabled 

For the direct CP asymmetry measurement we assume that the efficiency corrections and un- 
certainties cancel out. To account for uncertainty due to peaking background we use the following 
procedure. We assume an a priori flat distribution for the CP asymmetry between —1 and 1, which 
has a root mean square width of 0.58. We multiply this by the expected fractional contamination 
in our sample to obtain the systematic uncertainty. For (f)K^ {-k^ / rf) we assign 1.8% uncer- 

tainty, while for B~ K~ K~^K~^ we assign 3.5% uncertainty. For resonant B — > (j)K*{-^ Ktt^) 
events, the previous BABAR measurement JHl shows that the CP asymmetry is consistent with zero 
to within 9%. We therefore consider this to be negligible in our case. Using the control sample as 
we did with the branching fraction measurement, we vary the fixed input parameters of the fit to 
determine the uncertainty on the signal CP asymmetry. This was found to be 2.2%, bringing the 
total systematic uncertainty to 4.8%. 

In summary, we have performed the first BABAR studies of i? — > (j)K^ decay modes. The 
B^ (j)K^j branching fraction was measured, and an upper limit for the B^ (f)K^^ mode was 
determined. Our measurements are consistent with the assumption of isospin symmetry at the 
1.1% level. We have made the first measurement of the direct CP asymmetry in B~ — > <j)K~'j. 
For comparison, we quote the Belle results [7]: B{B~ (pK~j) = (3.4 it 0.9 it 0.4) x 10^^ and 
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1 „i 1 T^TTTT 1 

B^ (pK'^j 


5" 


Dource 


T T J. • J. 

Uncertainty 


T T J. ' J. 

Uncertainty 


K K^K 7 subtraction 


19.7% 


5.2% 


Feakmg Background 


6.4 % 


3.4% 


lit Bias 


2.6 % 


0.6% 


lit FDl parameters 


+7.0 0/ 
-5.9 /o 


+5.9 0/ 
-5.2 /O 


Total yield uncertainty 


11 u 

events 


17 ii 

-6 9 events 


Kaon Tracking 


2.8% 


4.2% 


Efficiency 


1.5% 


0% 


(j) Efficiency 


1.7% 


1.7% 


Particle ID 


2.8% 


4.2% 


hmgic Fiioton Efficiency 


1.8% 


1.8% 


Photon Spectrum Model 


U.4/o 


z.D/o 


Lo/Ln Cut 


1.2% 


1.2% 


7r°/?7 Veto 


1% 


1% 


Efficiency /acceptance uncertainty 


5.2% 


7.1% 


BB Counting 


1.1% 


1.1% 


Total 


+23 (V 
-22 /o 


+11.2(v 
-10.8 /o 



Table 2: Summary of the systematic uncertainties. 



B{B^ (pK^'-f) < 8.3 X 10-6 at 90% confidence level. The statistical uncertainties will improve 
as the B Factories collect more data over the next several years. In future measurements a large 
amount of CP violation would be a sign of physics beyond the Standard Model. 
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